Inelastic Interaction Corrections and Universal Relations for Full Counting Statistics 
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We analyze in detail the interaction correction to Full Counting Statistics (FCS) of electron 
transfer in a quantum contact originating from the electromagnetic environment surrounding the 
contact. The correction can be presented as a sum of two terms, corresponding to elastic/inelastic 
electron transfer. Here we primarily focus on the inelastic correction. 

For our analysis, it is important to understand more general — universal — relations imposed on 
FCS only by quantum mechanics and statistics with no regard for a concrete realization of a contact. 
So we derive and analyze these relations. We reveal that for FCS the universal relations can be 
presented in a form of detailed balance. We also present several useful formulas for the cumulants. 

To facilitate the experimental observation of the effect, we evaluate cumulants of FCS at finite 
voltage and temperature. Several analytical results obtained are supplemented by numerical calcu- 
lations for the first three cumulants at various transmission eigenvalues. 

PACS numbers: 73.23.-b, 72.10.-d, 72.70.+m 



I. INTRODUCTION 

Within the last years technological advancements have 
enabled the fabrication of sufficiently small (nanometer) 
solid state structures where electrons traverse the system 
as coherent quantum waves. The electron transport in 
such a quantum contact can be described with the scat- 
tering approach pioneered by Landauer and Biittikeri^. 
Within this approach, the contact is completely char- 
acterized by the set of transmission eigenvalues {T n }, 
< T n < 1, regardless its concrete structure. The (dif- 
ferential) conductance is given by the Landauer-Buttiker 
formula G = G Q E„ T n> G Q = 2e 2 /h. The same 
transmission eigenvalues determine superconducting and 
noise properties of the structure*^. Break junction 
experiments— provide excellent examples of tuning and 
experimental characterization of T„'s in concrete quan- 
tum contacts. In the scattering approach, the electron- 
electron interaction inside the contact is commonly dis- 
regarded. There is a good reason for that, eventually the 
same as for electrons in bulk metallic solids. Close to the 
Fermi energy, the only effect of interaction is to make a 
quasi-particle from an electron. These quasi-particles do 
not interact. This means that any contact at sufficiently 
low energies can be described within a non-interacting 
scattering approach. 

This however presumes an ideal voltage bias of the 
contact: the electrons are injected and absorbed by reser- 
voirs kept at a certain voltage. This assumption is too 
ideal: in reality, the contact is embedded in a macro- 
scopic electric circuit, and this electromagnetic environ- 
ment produces voltage fluctuations on the contact. The 
electrons traversing the contact can emit/absorb energy 
to/from the environment and interact by means of ex- 
change of photons that are present in the environment. 
The interaction due to the environment can not be dis- 
regarded at low energies and therefore becomes the most 



important interaction at low temperature and voltage. 
The environment is completely characterized by a fre- 
quency dependent impedance , in series with the quan- 
tum contact. 

The influence of the environment on electron transport 
has been studied in detail for tunnel junctions where all 
T n <C l 7 . The tunnel rates in the presence of an environ- 
ment can be evaluated for arbitrary impedance. For suffi- 
ciently large environmental impedances ZGq/2 ez> 1 
the interaction effects are large and the tunnel rates 
are strongly suppressed below a certain energy. This 
is termed Coulomb blockade of tunneling in a single 
junction 8 . The opposite case of small impedance z -C 1 is 
more realistic. In this case, the environment provides an 
interaction correction ~ z to the rates. This correction 
can be experimentally identified from its specific voltage 
and temperature dependence: it is seen as so-called zero- 
bias anomalyt^iS 

At arbitrary transmissions, the influence of the envi- 
ronment is more complicated and one cannot evaluate it 
for an arbitrary impedance (a progress in this direction 
has been reported inii). Still one can investigate the in- 
teraction correction ~ z to the contact conductance, both 
theoretically and experimentally. It has been demon- 
stratediSii^ that this correction is related to the second 
moment of current fluctuations: noise. The correction 
is proportional to shot noise ~ T„(l — T„), and disap- 
pears at perfect transmissions T n = 1. This prediction 
has been experimentally confirmed^* 

The environment influences not only the average cur- 
rent, but the whole statistics of electron transfers in the 
contact, the Full Counting Statistics (FCS). The theory 
of FCS for a quantum contact within the scattering for- 
malism has been developed malt. Later, it has been suc- 
cessfully applied to a variety of systems. The FCS is 
the statistics of current measurements over a given time 
interval r, r being much bigger than the typical time be- 
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tween electron transfers. It gives the probability P T (N) 
for N electrons to be transferred during this time inter- 
val. It is convenient to work with the generating function 
denned as 

F(V, X ) = J2 P -r( N ) eixN - (!) 

N 

The parameter \ is frequently called counting field since 
if one implements the Keldysh formalism for FCS^ x 
enters the formulation as a field conjugated to electric 
current. Derivatives of hiF(x) with respect to x gi ye an 
infinite set of cumulants of charge transferred, where the 
first two are related to average current and current noise. 

The environment influences the FCS in two distinct 
ways. A classical effect is governed by the impedance at 
low frequencies hcu <C max(eV, k^T) and scales like ZG, 
G being the total conductance of the contact. For the 
two first cumulants — average current and noise — this 
effect is nothing but voltage division between the contact 
and the external impedance. The effect is not that trivial 
for higher cumulants™, and the pioneering measurement 
of the third cumulan1ii& was shown to be affected by this 
environmental effect. However, this effect can be made 
arbitrarily small by a proper design of the low-frequency 
impedance to assure ZG <C 1. 

A more interesting effect comes from the impedance at 
frequencies Tiw > max(eV r , k^T) and scales as ZGq. This 
is the interaction correction discussed which eventually 
leads to Coulomb blockade effects at higher impedance. 
It has been shown ini^ that the correction can be sepa- 
rated into inelastic and elastic parts, the latter presents 
a renormalization of elastic scattering properties of the 
contact by interaction. It is feasible to observe the in- 
teraction correction to FCS in experiments, for instance, 
with well-characterizable break junctions. Such exper- 
iments would be certainly possible for the first three 
cumulants, and the developments in the fielci22i2ii22i2£ 
suggest that the higher cumulants can be accessed with 
proper measurement techniques as well. Apparently, the 
correction will include both the elastic and inelastic part. 
To provide theoretical support for these experiments is 
the main motivation and goal of the present work which 
concentrates on the inelastic part of the interaction cor- 
rection. 

In Section [H] we derive the interaction correction to 
the FCS in first order perturbation in z starting from 
the general form of the system-environment Keldysh ac- 
tion. The result allows for the identification of an elastic 
and inelastic contribution. In Section fllll we explain how 
one can re-interpret this result in terms of correlations 
of elementary events of charge transfer. We present con- 
crete analytical results in the next Section, particularly 
for vanishing temperature. 

While investigating the interaction correction in the 
limit of small voltages, we have found sets of simple re- 
lations for cumulants. Further analysis has shown that 
these relations are not specific for the setup considered 
and hold for any conductor regardless its properties and 



presence or absence of interaction. In fact, these univer- 
sal relations provide the generalization of the fluctuation- 
dissipation theorem^ for FCS. We discuss these relations 
at length in Section The derivation is provided in the 
Appendix. 

In the last Section we present our numerical results 
for arbitrary temperature and voltage. We have studied 
the experimentally interesting case of an RC-environment 
and plot the correction to the first three cumulants for 
conductors of different transmissions and a diffusive con- 
ductor versus eV/k^T as obtained from a numerical eval- 
uation of the integrals. In all cases we observe a crossover 
at eV ~ k^T which is related to a transition from ther- 
mal to shot noise behavior. 



II. ACTION 

The system we consider consists of a quantum conduc- 
tor which can be described by the set of its transmis- 
sion probabilities, T n , and a frequency dependent envi- 
ronmental impedance, z u in series. The voltage drop over 
the whole conductor-environment system is fixed. How- 
ever, the voltage in the node between the contact and the 
impedance can fluctuate. For instance, an electron trans- 
fered will momentarily charge the node creating a voltage 
pulse oc z that persists for some time and may influence 
further electron transfer. Thus there will be a fluctuating 
voltage in the node depending both on the probabilistic 
nature of the electron transfer in the quantum conduc- 
tor as well as on the impedance of the environment. We 
study the corrections to electron transport due to these 
fluctuations. 

It is convenient to work with the phase rather than 
voltage, which is defined as <f> = J dteV(t)/h. Since we 
study a quantum mechanical system, we have to describe 
quantum fluctuations of this quantity. This is most con- 
veniently presented in the language of Keldysh action 
that expresses physical quantities in the form of a path 
integral over the fluctuating phase on the Keldysh con- 
tour. Since the Keldysh contour consists of two parts, 
the integration proceeds over two sets of variables (j)^ (t) 
corresponding to these parts. 

The Keldysh action approach to mesoscopic quantum 
circuits has been pioneered in2£ and has been extended 
to cover FCS and arbitrary quantum contacts i n 17 i 26 . As 
usual in the theory of FCS, it is the generating function 
F(x) which is presented as a path integral over the fluc- 
tuating phase. The action in the path integral is a sum of 
the actions describing the constituent parts of the circuit: 
the conductor action, S c , and the environment action, 

F( X ) = /'# + ^-exp{-iS'[^ + ,^-]}; 

S[<j> + ,^] = S c {cf > + ,4>-} + 

+ S cnv [4> + X /2-0 + ,<I>-x/2-<r]. (2) 
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We use superscripts ± to denote the phases at dif- 
ferent parts of the contour and use traditional no- 
tations if (x) f° r their half-sum (difference). Cur- 
rent, noise and higher moments of FCS follow as 
dhiF/dx\x=Q) <9 2 ^ n F/d X 2 \ x =o, ' • '■ For a linear environ- 
ment the action S onv is a bilinear function of the phases 
depending on the impedance and temperature T only, 



S _J_ 

-(t),cf>-(t))A(t-t') 



dt / dt' 



0+(f) 
<t>-{t>) 



(3) 



where the coefficient matrix, A, depends solely on fre- 
quency and temperature: 



A{w) = 



-iu[z- x + 2NJ&ZZ 1 ] 
~[A+-(-u)]* 



-[A++(-u)Y 



(4) 



N u = {exp[huj / ksT] — is the Bose-Einstein distri- 
bution function. 

The most concise way to write the conductor action is 
in terms of Keldysh-Green functions L,Roi the left /right 
reservoii»i&: 



l + ^({L,R}-2) 



(5) 



The trace is over Keldysh space and energy/frequencies. 
In equation (J5J the fluctuating phase in the node enters 
in the form of a gauge transform in one of the reservoirs 
as L = e^G rcs e-^, R = G rcs , where = tp + x/2r z . The 
equilibrium Keldysh-Green function depends on time dif- 
ference (or energy) only and reads in a usual way 



G 



1 - 2/(e) 2/(e) 
2(1 -/(e) 2/(e)-l 



(6) 



with / being the Fermi distribution function of the cor- 
responding load. Equations 12I3I4I5I6H define our model 
and all we have to calculate is the cumulant generating 
function F for a given environment z w and conductor T n . 
However in the general case this is a formidable task. A 
natural way to proceed is to assume z <C 1 and thus 
to treat the effect of the environment as a perturbation. 
In zeroth order (no environment), the phases are obvi- 
ously related to the applied voltage and do not fluctuate: 
<p = eVt. 

Putting this in equation J5| and taking the trace gives 
the well known expression for the generating function 
of a mesoscopic conductor in terms of its transmission 
eigenvalues^ 

n J 

ln{l+T n (e^ - - f r ) + {e- lx - - fx)} (7) 

Here and in the following indices 'l,r' refer to the 
left /right lead. We assign the voltage to the left lead, 
so that fi(e + V)=f r (e) = f(e). 



The first order correction is proportional to the fluctua- 
tions of 4>, which are small oc z. Expanding the logarithm 
in equation (JjJJ to second order in <fi gives a second order 
contribution 

-iSSW = Zpn{ZM(2)}_I (^j\v{DA^DA^) 

= ^DTr{A ( V}-^ (J^j DD+Tr{A^A^} 

where the following relations hold under the trace 

= A^ 2 - A4>L + (j)L - 2cj)L + R + (j)LR 

—2(j)R + L + (j>RL + 4(j)L + R + L + (t)R, (10) 
(j)L + cj)R + cj)R + (t)L - cj) 2 (RL + LR) (11) 



A$ 



and 



D- 1 = l + T n [(e^-l)/ i (l-/ r ) + 
+( e -«X-l)/ r (l _/,)]. 



(12) 



All quantities with superscript '+' are taken at energy 
e + oj and integration over energy and frequency is im- 
plied. For convenience we omitted the explicit depen- 
dence on £, u>. It is easily found from the definition of the 
trace. The first term in equation i|ll|) for instance would 
read T„/8Tr J d£dujD £ (j>^ U jL e+U j(j> UJ R e , where the trace is 
understood over Keldysh indices. 

These terms are quadratic in phase, and by virtue of 
the path integral in equation are to be replaced with 
their averages given by the environmental action. These 
averages read 



m vx 
v*x \x\ 2 



(2N„ + 1) 







(13) 



After some ordering of terms, the first order correction 
to the cumulant generating function can be presented as 



(lnF( X ))M 



T 

ft JO 



du> 



(2N U + 1)S£ ) ( X )+ 



N u S£>(u,x) + (N u + l)s£>(-<".x)J ■ (14) 

We note that there are three different terms which can 
be identified as being due to elastic electron transfer, and 
inelastic transfer with either absorption (positve uj) or 
emission (negative u) of energy Tim from/to the environ- 
ment. Explicitly, in terms of filling factors these terms 
read 



5 el = -2^r n (i-r„) 



dT n 



(15) 



E 



J^-{DD+[T n (f l -f+) 

+2T n (e«-l)/,(l-/+) 
+4T„ 2 (cos X -l)//(l -/,+)(/+ 

+(1-D)(1-T„ -£>+)} + 1 



~fr)} 
I <-> r 



X 



-x 



(16) 
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Since the expression is symmetric with respect to ex- 
change I <-> r of the leads and simultaneous change of 
the sign of the counting field, the cumulants are either 
even or odd functions of the voltage applied. Following 19 
we present the elastic part of the correction as a change 
oc z of transmission eigenvalues. It is clear that the main 
contribution to the inelastic part of the integral in equa- 
tion (|14|l comes from small frequencies less than the max- 
imum of voltage and temperature. Elastic processes do 
not have this restriction and S- m contributes to the inte- 
gral at all frequencies. 



III. INTERPRETATION 

The advantage of the FCS approach to quantum trans- 
port is that in many cases the FCS expression can be 
re-interpreted in terms of elementary events thus provid- 
ing some insight into the relevant transport properties. 
A well-known example of such interpretation is provided 
by the non-interacting Levitov formulas^. At vanishing 
temperature, it allows to present the statistics as a super- 
position of reV/h elementary "games" in each transport 
channel n , each "game" resulting in either transmission 
(with chance T n ) or reflection (with chance 1 — T n ) of an 
electron coming to the contact. A recent example of such 
a re-interpretation concerns spin statistics 2 ^. In this Sec- 
tion, we show how to interpret the inelastic interaction 
correction given by Eq. 1161 As we will see, the interpre- 
tation is probably too complicated to be constructive. 
Still it gives some insight to the form of the expression, 
in particular, the presence of denominators D,D + . 

To give an interpretation of the form of S m let us 
go back to relation Q that holds for non-interacting 
electrons. In the limit of large measurement time r, 
Ti/t <C max(/cBT, eV), we can discretize the integration 
over energies. The general structure of the generating 
function 28 is then 



^ (o) (x)=n A ( i? )' 



(17) 



with X(E) = D~ l . The product is taken over dis- 
cretized energies as well as over transport channels. The 
right hand side of equation (|T7|l is a product, so each 
term in this product can be regarded as an independent 
process ("game"). The concrete form of X(E) suggests 
that three distinct outcomes of each process are possi- 
ble: (i) electrons transmitted from the left reservoir to 
the right with probability P + = T n /;(1 — f r ), (ii) trans- 
mission from right to left with P_ = T n f r (l — /;), (iii) 
no transmission (Pq = 1 — P + — P-). Indeed, the gen- 
erating function for each process is then A — ^ Q P a X a 
with X a — e lax . Electrons at different energy (and/or 
channel) arc uncorrelated since the complete generating 
function factorizes in terms of X(E). 

How can this picture change if one introduces electron- 
electron interaction via the environment? It is clear that 



interaction will bring about all kinds of correlations be- 
tween electrons at different energies and the simple pic- 
ture presented above is not true anymore. The major 
change is, that the factorization in uncorrelated elemen- 
tary events does not hold. Presumably, the generating 
function of an elementary event A will depend on many 
different energies. The relation (|16|) suggests that in low- 
est order in z it depends on two energies only, A(E, E'), 
where \E—E'\ = hui is the energy of an absorbed/emitted 
photon. With this accuracy, the cumulant generating 
function can be expressed as a product over pairs of en- 
ergies 



F (o) (x) + F (D (x)= TJ A(E, E'). 



E,E> 



Without interactions, 



A(E, E') 



X(E) if E = E' 
1 if E ^ E' 



(18) 



(19) 



so that electrons with different energy are uncorrelated. 
If 6A(E, E') is the interaction correction to A, the change 
of the cumulant generating function reads 



(ln^(x)) (1) = E 



E.E 



6A(E,E') 
, X(E)X{E') 



(20) 



This explains already the presence of denominators in 
()16J) . In addition, we conclude from (|16(l that contains a 
single sum over transport channels, that the elementary 
events do not involve electrons in different channels, even 
though they involve electrons at different energies. This 
is probably valid only for the first order correction. 

To proceed, let us note that the correction SA(E, E') 
consists of terms to be divided into three classes. Firstly, 
there will be terms presenting new-events, SA now (E, E') 
not taking place for non-interacting electrons. An exam- 
ple is an electron transfer from the left to the right with 
photon emission. As we can assert from (11611 it comes 
with probability oc T n fi(E)(l - f r (E'))(l + N E -e>). 
Another example is a two-particle process consisting of 
elastic electron transfer at energy E accompanied by 
inelastic transfer, its probability being proportional to 
T*ME)f r (E)(l - fi{E>)){l - f r {E)){\ + N E - E ,). Sec- 
ondly, since the probabilities of elementary old- events are 
modified by interaction, there will be terms depending on 
a single energy only, those can be seen as the modification 
of X(E), X(E) -> X(E) + 5A old {E,E). They are incorpo- 
rated into the elastic part of the correction. Finally, the 
environment will introduce correlations among pairs of 
old-events, represented by 8A COII (E, E'). For instance, 
the correlation between left-right transfer at energy E 
and right-left transfer at energy E' will come with a fac- 
tor Tlh(E)(l - f r (E))f r (E')(l - ME')). _ 

These three contributions simply add up in the correc- 
tion to the generating function, 



(lnF( X )) 



(i) _ 



E 



SA old (E,E) 
X(E) 



5 



; E 

E>E' 



5A ncw (E,E r ) + 5A cm (E,E>) 
\(E)\(E') 



(21) 
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One recognizes this structure in the relations for S c \ (first 
term) and for S- m (second term). In principle, in this way 
one can recover the correction to the generating function 
of an elementary event SA and find the (corrections to) 
probabilities of all possible outcomes, new ones as well as 
old ones. This gives a re-interpretation of the correction: 
any term of equation i|16[l is assigned to a term of one of 
the three classes discussed. 

However, the procedure is cumbersome and hardly 
practical because of the large number of possible pro- 
cesses and outcomes. For a two-electron process, each 
incoming electron can be in one of four possible states 
(coming from the left or the right at E or E'), the same 
for outgoing electrons. This gives in total 2 8 terms: it 
looks like a somewhat lengthy interpretation of a rela- 
tively compact equation l|ltjf> . This prevented us from 
accomplishing this program explicitly. We are satisfied 
with the fact that the combinations of electron and pho- 
ton filling factors make sense for the terms contributing 
to <5A ncwcorr . The picture and the interpretation are ex- 
pected to become even more involved for the corrections 
of higher orders in z. 



IV. ANALYTICAL RESULTS 

At vanishing temperature we can easily perform the 
integration over e in equation (|14|) . The full correction 
to FCS then reads 

eii t f°° Vz 

(lnF( X )) (1) = -eVS cl / dco^- (22) 

S ol = -2^T„(l-T„) r 



eV/h 

e ix - 1 



+ T n (e < x-1) 



It has a simple and interesting structure revealing the 
relationship between the elastic and inelastic part of the 
correction. If we took the elastic part only, by virtue 
of J7| we would obtain a similar expression. The differ- 
ence is that the integration over u) would start at zero. 
Therefore, the inelastic part of the correction precisely 
cancels the modification of the elastic transmission for 
photon energies in the interval < hto < eV. It is indeed 
expected from general reasoning^ that the low-energy 
divergences present in the elastic correction are cut off 
at energies ~ eV. Our somewhat unexpected result is 
that at vanishing temperatures this cut-off is sharp and 
clear. In agreement with expectations, the inelastic S{ n 
only contributes at frequencies < Tluj < eV reflecting 
the fact that the only energy source for inelastic processes 
is given by the voltage. 

The correction to the rn-th cumulant S^ m ^ is given by 
derivatives of the above relation, we find: 



ss (m) 



T 

1 n. 



d(i X ) m l + T n (e^-l) 



JeV/h ^ 



Importantly, the correction is proportional to the ro+l-th 
cumulant for the non-interacting case, 
This generalizes the resultsiSiii for the average current. 
The environment enters the corrections as an integral 
over the impedance and affects every cumulant in the 
same way. 

A common and experimentally interesting model for 
a frequency-dependent impedance is that of an RC- 
environment, z u = z(l + iw/o; c ) _1 . The impedance is 
cut at uj c — 1/ (RC) and approaches a constant value 
of z at 10 -C lj c . The integral governing the correction 
evaluates to 



duj^^- = In \ 1 



eV/h 



e 2 V 2 



Tiuj c 

In , if huj c 3> eV 

eV 

That is, it diverges logarithmically at sufficiently low 
voltages eV -C %u) c . This is the well known zero bias 
anomaly. As has been shown it holds for any cumulant. 
Different cumulants have the same functional dependence 
on voltage and can be scaled by the prefactor of equation 
(|24|l which depends only on the transmission probabili- 
ties. 

Even if z <C 1 the correction oc zhi(huj c /eV) can be- 
come big at sufficiently small voltages, eV/hu; c ~ e~ 1/>2z . 
It has been shown in — that in this case one has to im- 
plement the renormalization procedure neglecting the in- 
elastic part. The elastic correction can be consolidated 
in the energy dependent renormalization of transmission 
eigenvalues given by 



dT n (E) 
d\ogE 



2zT n (E)(l-T n (E)). 



(25) 



We do not consider this here. Rather, we expect that 
finite temperature will lead to a rounding off of the singu- 
larity at small voltages in the same way as for the current 
correction. 

The equation l|14|) is too complex at finite tempera- 
tures so that it is hard and non-instructive to perform 
the integration over energies e. However, the correction 
to any cumulant of finite order derived from the gener- 
ating function is an integral over a finite polynomial of 
Fermi functions. This integration can easily be done for 
arbitrary temperature and voltage. The analytical for- 
mulas obtained in this way are too long to be of any use 
except numerical evaluation. The correction to any cu- 
mulant is proportional to T n (l — T n ), thereby vanishing 
at perfect and vanishing transmission. 

For the correction to the average current we find 



51 



-2eJ2 T n{l-T n ) 



du> 



^ h ^T-f sinh^ 



CO 



cosh -rXr 



cosh ^ 

Kb J 



(26) 
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I All; 



(27) 



where (|A7|^) is the finite-frequency current correlator 
without environmenti^S. The correction to the current 
is thus related to the noise in the absence of an environ- 
ment. 

The correction to the m-th cumulant is an m + 1-th 
order polynomial in T n . The term linear in T n has the 
same functional dependence as that in (|27l) . The reason 
for this is that in the limit of small T n the Full Counting 
Statistics is that of a tunnel junction: electron transfers 
are rare and consequently independent. We get a super- 
position of two Poissonian statistics for electrons tunnel- 
ing from the left to the right and from the right to the 
left expressed as 

lnF(x) = t (r LR (V,T)(e l * - 1) + T RL (V, T)(e"« - 1)) 

^lr,rl being tunneling rates in these two directions. 
The interaction in this limit modifies Tlr^rl, this be- 
ing the only effect on FCS. We will see below that these 
two rates are related by the detailed balance condition 
Trl = Tlr exp(— eV/k B T). From this it follows that in 
the tunneling limit 



_g{m) _ 
T 



/coth(^) 



if m odd 
if m even 



(29) 



for any interactions. 

Analytical work in the limit of small voltages eV <C 
k B T gave us some relations between the cumulants. How- 
ever, we have recognized that these relations are not spe- 
cific for the interaction correction but are of general na- 
ture. That is why we discuss them in the next Section. 



V. UNIVERSAL RELATIONS FOR 
CUMULANTS 

The detailed investigation of the interaction correc- 
tions to FCS is hardly effective without appreciation of 
universal relations for FCS cumulants that hold with no 
regard for interaction and/or concrete structure of the 
conductor. This is why in the course of this research we 
had to understand the general constraints imposed on 
the FCS by laws of quantum mechanics and thermody- 
namics. 

We show in the Appendix that this universal relation 
can be most generally expressed in the following form 



F(V,x) = F(V,- X + ieV/k B T). 



(30) 



A didactic representation of this relation can be obtained 
by recalling the definition of F as generating function of 
the probability distribution of a certain number N of 
particle transfers, 



F(V, X )^J2 P ^ XN - 

N 



(31) 



Applying i|30[) . we observe that the probabilities of oppo- 
site number of particles transferred are related by 



Pn(V) = P_ N {V)e eVN / k * T . 



(32) 



This relation is well-known for independent tunnelling 
events (see, e.gJ-^S.) and is referred to as detailed balance 
condition. We thus demonstrate that the detailed bal- 
ance holds for any N irrespective of possible interactions 
and correlations in and beyond the conductor. 

Whatever didactic, the detailed balance condition is 
not easy to apply to cumulants. We do this with the 
universal relation l|3U|) that obviously holds for InF as 
well. A series of relations for cumulants is obtained by 
taking derivatives of this relation with respect to voltage 
at V = 0. First of all, we just go to the equilibrium limit 
V — > to obtain 



lnF(x) =lnF(- X ), 



(33) 



all even cumulants thus vanish at equilibrium. This rela- 
tion and all subsequent relations till the end of this Sec- 
tion are valid only in the limit V — *■ 0. 

Taking the first derivative with respect to voltage, we 
arrive at 

9 (lnF( X )-l n F(- X )) = -^^M. (34) 



deV 



If we expand this in x, we obtain a series of equations 
that relate even cumulants with voltage derivatives of 
odd cumulants 



2k B T dS( 2n+ V 
~e dV 



(35) 



The first equation in this series is nothing but Johnson's 
noise formula, 



e 2 ^ = 2fc R T 



(36) 



that relates zero-voltage conductance and equilibrium 
current noise. 

The next series is obtained by taking the second deriva- 
tive with respect to voltage and making use of (1341) . 



d 2 
dV 2 
e 8 2 



k B T dVd X 



(lnF(x)-lnF(- X )) = 
(lnF(x) + lnF(-x)) 



(37) 



which is only practical if even(odd) cumulants are not 
even(odd) functions of voltage, that is, in the absence of 
electron-hole symmetry. Since our model is electron-hole 
symmetric, this relation is of no immediate relevance. 
The first relation in the series reads 



d 2 i 

W 2 



ds^ 2 



k B T dV 



(38) 



It relates d.c. current induced by low-frequency a.c. volt- 
age (rectification effect) to low-frequency current noise 
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proportional to d.c. voltage applied. This relation was 
discussed in detail in^fi in the context of the photovoltaic 
effect. 

Taking the third derivative with respect to voltage and 
making use of 1(38(1 we obtain another series: 

53 n jrf \ 1 j?t w ie3 d 3 lnF(x) 
W 3 (lnF(x) lnF ^ = (fc B T)3 9 X 3 



2>ie d 3 



(lnF(x)+lnF(- X )) (39) 



k B T dV 2 d X 

The first relation in the series can be rewritten as 



2e 



,d 2 SW 
dV 2 



1 

3 V dV 3 



(40) 



{k B T) 2 

The l.h.s. gives the change of the current noise induced 
by low-frequency a.c. voltage at non-matching frequency. 
This response coefficient, and its importance, has been 
recently discussed in£i, where it has been termed "noise 
thermal impedance" . The authors have conjectured the 
relation of this coefficient to the fourth cumulant. 

It is easy to see that Eq. (|30|) holds for the elastic part 
of the FCS even before integration over energies in J7J, 
since for any energy 



flQ f r ) _ „eV/k B T 



(41) 



/r(l"/j) 

The corresponding proof for the inelastic contribution 
can be done only after integration over energies, is cum- 
bersome and has provided a good check for the validity 
of expression 116(1 . 

VI. NUMERICAL RESULTS 

We restrict the numerical analysis to the RC- 
environment model discussed above. There are three en- 
ergy scales in the system: voltage, temperature and the 
cut-off frequency of the environment, u> c . It is especially 
interesting to study the case huj c 3> eV, k^T. In this 
case it is expected that all cumulants show the logarith- 
mic divergence at small voltages/temperature, which is 
well known for the conductance: SG ~ Go In rrh t\ ■ 

^ max(e V,k&T) 

This motivated us to the following choice of the presen- 
tation of the results: we plot the second derivative of 
the correction to the first three cumulants with respect 
to voltage versus eV/k^T ( figures Hl2l3|) . It is the sec- 
ond derivative that approaches a limit independent of uj c 
upon increasing the ratios hui c /eV, Titu c /kB,T. To illus- 
trate the dependence on transmission and to assess differ- 
ences between specific conductors (ballistic, tunnel,. . . ), 
we plot the results for a single-channel conductor with 
transmission values ranging from 0.1 to 0.9. Another in- 
teresting reference system is a diffusive conductor. The 
results for a diffusive conductor can be obtained by av- 
eraging over transmission eigenvalues with the distribu- 
tion function p(T n ) = (T n -\/l — T n ) Since all cumu- 
lants are polynomials in T n this is equivalent to replacing 
T n^ v^T(™)/r(m + l/2). 



> 

(D 



C/3 

CM 




eV/k B T 



FIG. 1: Second derivative with respect to voltage of 
the correction to the average current, d 2 S {1) /d(eV) 2 = 
hk B T/[2eRG Q J2 n Tn(l - T n )]d 2 1 /d(eV) 2 vs eV/kaT. 



At eV/k-oT 3> 1 the corrections are dominated by S e i 
and the emission term (— cj, x)i since the environment 
does not provide any energy at T — > 0. The functions 
plotted oc z/V for all cumulants, and are given by Eq. 
123 As discussed in the paragraph below that equation, 
this leads to a suppression of the conductance (and any 
other cumulant) at small voltages which is termed zero- 
bias anomaly. 

Let us now discuss differences between cumulants, 
starting with the current (fig. An apparent feature of 
the current correction is that the corrections to different 
conductors can be scaled to one curve by the common 
prefactor ^2 n T n (l — T n ), (Eq. |2EJ. The corresponding 
curve for a diffusive conductor can be obtained, follow- 
ing from averaging over transmissions, by multiplication 
with 2/3 (and removing the dependence on T n in the 
normalization). This feature is unique for the correction 
to the current and independent of the choice of a spe- 
cific environment. Higher order cumulants have a more 
involved dependence on transmission eigenvalues and en- 
vironment. 

In the opposite limit of large temperatures, eV <C k^T, 
the environment provides the energy for inelastic electron 

transfer. Consequently the absorption term, S^\u>,x) 
becomes more important in equation 1141 It is expected 
that in the same way as for the zero-bias anomaly the 
correction to the conductance is logarithmically diverg- 
ing with temperature. Due to the choice of presentation, 
this term is not visible in figure ^ What is shown in 
this plot and the following is the lowest order term in the 
expansion in powers of eV/k^T. For the current (and 
any odd cumulant) this is a linear term, following from 
the symmetry with respect to inversion of voltage as ex- 
plained below equation 1161 In summary, the correction 
to the current shows a crossover from a temperature to 
voltage dominated regime at eV ~ k^T. 

In figure [21 we plotted the second derivative with re- 
spect to voltage of the correction to the noise for several 
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single channel conductors and a diffusive conductor. At 
small voltages, all curves start with zero slope since noise 
is an even function of voltage. Interestingly with increas- 
ing voltage they all cross the x-axis at eV ~ k^T, before 
approaching zero. As expected, changes occur on the 
scale of temperature. The limit of vanishing tempera- 
ture, eV/k^T 3> 1, namely the proportionality to z/V, 
can be discussed along the same lines as for the current. 
However there are several striking differences that were 
absent in the correction of the current. 

Unlike for the current, curves for different transmis- 
sion can not be reconciled by scaling, rather we observe 
a strong dependence on the value of the transmission. 
This dependence is non-monotonous. However at small 
voltages, conductors with T n < 1 have positive correc- 
tion while those with T n <C 1 have a negative. This 
sign-change could have been conjectured since it is well 
known that corrections to the m-th cumulant are related 
to the unperturbed m + 1-th cumulant. Hence, the cor- 
rection to noise should be related to the third cumu- 
lant whose dependence on the transmission eigenvalues 
(oc T„(l — T n )(l — 2T„)) changes sign at intermediate 
transmissions. 

This behavior can be interpreted by looking at the ex- 
treme cases of T n <C 1 and T n < 1. It is plausible that 
either conductor in the absence of an environment pro- 
duces little shot noise since in the first case the current 
is "most of the time" zero with only rare transfers of 
charges. In the second limit, electrons are transfered with 
probability close to one and only occasionally reflected. 

The same conductors embedded in an environment 
however will feel a suppression of current due to fluctua- 
tions of the voltage in the node as discussed in sectionlTTl 
For the noise of the tunnel conductor that means, that 
the rare electron transfers being the source of noise, will 
be suppressed furthermore leading to a reduction of noise 
(negative correction). If T n < 1, the suppressed conduc- 
tance means that the quasi-constant flow of electrons will 
be interrupted more often, resulting in an enhancement 
of noise (positive correction). Consequently there will 
be a crossover at intermediate values of the transmission 
which can be seen in figure ^ 

Comparing the shape of the curves we observe that 
the maximum (absolute) value lies for most conductors 
at V = 0. Remarkable exceptions are the diffusive con- 
ductor and such with transmission close to the crossover 
(T„ = 0.4, inset figure We might argue that the dif- 
fusive conductor inherits this feature from intermediate 
transmissions or rather that it is an "interference" ef- 
fect determined by the coefficients of lowest and highest 
power in transmission (T n ,T^) in the expression for the 
correction to the noise. 

The corresponding plot of the correction to the third 
cumulant, which reflects the asymmetry of electron trans- 
fer, is presented in figure^ It shares features of both cur- 
rent and noise correction. Due to the different symmetry 
with respect to voltage inversion, the corrections to the 
third cumulant start at zero with linear slope. Again 
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FIG. 2: Second derivative with respect to volt- 
age of the correction to the noise, d 2 / d(eV) 2 = 
hk B T/[2e 2 RG Q ]d 2 I (2) /d{eV) 2 vs eV/k B T. The general 
tendency as a function of transmission T n is indicated by an 
arrow. The inset shows a zoom for T n = 0.4. 



there is a crossover at eV ~ k^T and a decay with z/V 
at large voltages reflecting the zero bias anomaly. The 
dependence on T n is non-trivial which is not surprising 
since the expression contains four terms of different power 
in T n , each of which can have a distinct dependence on 
eV/k-QT. However ballistic/tunnel conductors at small 
voltages separate in the same way (albeit with opposite 
sign) as for the correction to noise to the lower/upper 
part of the plot. Comparing curves of small T n with those 
for the current we recover the tunneling limit discussed in 
section llVl Interestingly the correction for intermediate 
transmission appears to be more feature-rich than for the 
limiting cases. This is clearly an indication that terms of 
different power in T n contribute with equal weight to the 
interaction correction while the correction for a tunnel 
contact is dominated by just one term. 




FIG. 3: Second derivative with respect to voltage of 
the correction to the third cumulant, d 2 S^ /d(eV) 2 = 
hk B T/[2e 3 RG Q ]d 2 I {3) /d(eV) 2 vs eV/k B T. The general ten- 
dency as a function of transmission T n is indicated by an 
arrow. The inset shows a zoom for T„ = 0.6. 

As the main result of our numerical analysis we note, 
that the corrections to cumulants strongly depend on the 
transmission of the contact. They can have either sign 
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and a distinct dependence on the ratio of voltage and 
temperature. Both of these facilitate the experimental 
detection of environmental effects on transport proper- 
ties. The plots in this section were obtained for an RC- 
environment. In principle one could obtain results for 
any given z w . At least qualitatively we expect the cor- 
rections due to a (physical) environment to be similar to 
those presented. 



VII. CONCLUSION 

We have studied the interaction correction to Full 
Counting Statistics of electron transport in a quantum 
contact. It was shown that the interaction can be mod- 
eled by an environmental impedance in series with 
the contact. In section [H] we presented a formulation of 
the problem in terms of a non-equilibrium Keldysh ac- 
tion. Assuming ZGq C 1 we proceeded perturbatively 
and calculated the correction to the cumulant generat- 
ing function, that is, to any cumulant (Eq. I14|) . This 
correction splits into three parts corresponding to elastic 
electron transfer and inelastic transitions with absorp- 
tion/emission of energy from the environment. 

We looked in detail at the structure of the interaction 
correction and found a re-interpretation in terms of ele- 
mentary events. This provided a deeper insight into the 
physics involved and presented a basic check for the ob- 
tained expression. Since the full expression, Eq. 1141 is a 
complicated function of temperature and voltage which 
is not easily understood, we looked at certain limiting 
cases. In the limit of vanishing temperature we found a 
particularly simple expression for the correction to any 
cumulant, Eq. 1231 For the opposite limit of vanishing 
voltage, we realized that any expression between cumu- 
lants is due to a universal relation of detailed balance 
for the generating function that holds irrespective of the 
concrete structure of the quantum contact and possible 
interactions, Eq. [3U] In order to bridge those limits and 
to enable the experimental observation of environmen- 
tal effects on electron transport in a quantum contact, 
we calculated numerically the correction due to an envi- 
ronment to the first three cumulants for arbitrary volt- 
age, temperature and different transmission eigenvalues. 
We have shown that the corrections show an interest- 
ing crossover behavior from voltage to thermal noise at 
eV ~ IcbT as well as a specific non-trivial dependence on 
transmission eigenvalues. The presented analytical and 
numerical results facilitate the measurement of the inter- 
action correction. 
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APPENDIX A: DERIVATION OF THE 
UNIVERSAL RELATION FOR FCS 

In this appendix, we will present the derivation of 
the universal relation (|3L)[) that is valid for FCS of 
any conductor regardless its concrete realization. A 
well-known example of a relation of this type is pro- 
vided by the fluctuation-dissipation theorem that re- 
lates the linear response and Gaussian fluctuations. 24 
The same approach can be extended to provide similar 
relations for non-linear frequency-dependent responses 
and non-Gaussian fluctuations22ii2ii2i using Hermicity, 
time-reversibility and KMS^2& relations. These results 
cannot be immediately used for our purpose since they 
are formulated in terms of relations between multi-point 
Keldysh Green functions rather than in terms of gener- 
ating functions. 

We apply this approach to the most general generating 
functional of current fluctuations where both voltage ap- 
plied V and counting field x depend on time. The deriva- 
tives of the functional with respect to x(t) gi ye averages 
of current operators. The Hamiltonian in the presence of 
the voltage source can be written as 

H(t) =H - m{t)I/e; = eV(t)/h, (Al) 

/ being the operator of full current in the conductor. 

We make use of the interaction picture introduc- 
ing I(t) = e iH t/hj e -iH t/h^ The genera ting function 
reada 16 . 17 ! 26 

F({*+(t)},{0-(t)}) = <L>t( {rW} )f7( W +(t)})),(A2) 
U(W + (t)}) = fexp (i J dt <t>+[t)i{t)/^j , (A3) 

UH{<t>-(t)}) = ^exp (-i J dt 0-(t)J(t)/e) , (A4) 
(...)= Tr (. . . p) ■ p = e-"/ k * T /Tr (V^ bT ) ,(A5) 

where ~T(^T) stands for (anti)time-ordering of the oper- 
ators and ± (t) = (p{t) ± x(*)/2- This expression is for- 
mally equivalent to the generating functional for multi- 
point Keldysh Green functions used in 33 '24. The only dif- 
ference is that the Green functions generated are those 
of current operators. 

We shall assume time-reversibility of the Hamiltonian. 
Since in this case 

H T = H; I T = -I; F{t) = -I(-t) (A6) 

we observe the following transposing rule for S'-operators: 

(u(mm) T = #({-#-t)» (a?) 

Transposing operators in the average (|A2|1 . we obtain 

F({4>+(t)},{<j>-m) = (u({-0 + (-t)})u\{-r(-tm). 
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In comparison with (|X2j, U,W are interchanged. We 
want them back to their original positions. The way to 
do this is to make use of KMS relations: For any operator 
A 

pA{t) = A(t + iti/k B T)p. (A8) 

We do this commutation with all operators comprising 
to obtain 

pua-^i-t)}) = 

= T*exp(-ifdt<l) + (-t)I(t + ih/kBT)/eJp . (A9) 

We shift now the time argument of <p + by ifi/k B T to 
obtain 

"fexp (-i J dt<f> + {-t)i{t + ih/k B T)/e^J = 
U({-</>+(-t + ih/k B T)}). 



This step looks rather heuristic. Since nothing is assumed 
concerning the analytical properties of as a function 
of complex time, the complex shift may be ambiguous. 
However we note that we are mainly interested in <f^ (t) 
that change at time scales much bigger than h/k B T: for 
those, we expect no ambiguity. 

Finally, we cycle operators under the sign of trace to 
obtain 

F({0+(t)},{^-(t)}) = 
F({-<p+(-t + i%/k B T)}, {-4>- (-*)}). (A10) 

For quasi-stationary V, x we substitute <j> = eVt/K ± 
X/2 and neglect the dependence on time-independent 
phase to arrive to 

F(V, X ) = F(V,- X + ieV/k B T) (All) 

which is the universal relation to prove. 
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